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Abstract. We consider open systems generated from one-dimensional maps 
that admit a finite Markov partition and use the recently developed theory of 
isospectral graph transformations to estimate a system's survival probabilities. 
We show that these estimates are better than those obtained through a more 
direct approach. 

1. Introduction 

Recently, a nontrivial relation between the dynamics of networks and the dy- 
namics of open systems has been found. This discovery has lead to advances in 
the analysis of network dynamics and has also introduced a new research direction 
concerned with the finite time properties of open systems with finite sized holes. 

The main idea behind this approach is that network dynamics can be broken 
down into three parts; (i) the network's graph structure, often referred to as its 
topology, (ii) the local or intrinsic dynamics of the network elements, and (iii) the 
network interactions between these elements. To each of (i)-(iii) there is an associ- 
ated dynamical system, which together can be used to characterize the dynamics 
of the network [21 El O [B] ■ The same approach can be used to study the topological 
properties of open dynamical systems [TJ || • 

One of the fundamental concerns in the study of networks is understanding the 
relation between a network's structure and its dynamics. However, the networks 
we often encounter in either nature or engineering are typically very large, i.e. have 
a large number of elements. It is therefore tempting to want to reduce such net- 
works by excluding some subset of these elements while preserving some important 
characteristic(s) of the original network. 

Since real networks are typically dynamic, one might first consider the spectrum 
of the network's weighted adjacency matrix to be one such a characteristic worth 
preserving. That is, one could hope to find a way of reducing a network while main- 
taining its spectrum. However, it seems immediately clear that such an "isospectral 
reduction" is impossible based on the Fundamental Theorem of Algebra. 
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In fact, it is possible to reduce the size of a network (or matrix) while preserving 
its spectrum. This theory of isospectral reductions can be found in [B] and is 
part of a larger theory of isospectral transformations introduced in the paper. Such 
transformations were used in [7j to improve each of the classical eigenvalue estimates 
associated with Gershgorin, et. al. [ID] . 

In the present paper we make use of this interplay between dynamical networks 
and open dynamical systems. We obtain estimates of survival probabilities for 
a class of one dimensional maps which admit a finite Markov partition. Most 
importantly, we shown that any isospectral transformation corresponding to an 
open system leads to sharper estimates of system's survival probabilities. 



Let / : I — > I where I — [0, 1]. For = go < <Zi < ■ • ■ < Qm-x < <lm = 1 we 
let £j = for 1 < i < m and assume that the following hold. First, the 

function is differentiable for each 1 < i < m. Second, the sets & = (<fo-i,<&] 
form a Markov partition £ = {^ij^i of /. That is, for each 1 < i < m the closure 
cl(f(£i)) is the interval [qj,qj + k] for some k > 1 and j that depends on i. 

We consider the situation where orbits of / : I — s- I escape through an element 
of the Markov partition £ = or, more generally, some union H of these 

partition elements. Equivalently, we can modify the function / so that orbits cannot 
leave the set H once they have entered it. Here, orbits that enter H are considered 
to have escaped from the system. This later approach of modifying f\u turns out 
to be more convenient for our discussion and will be the direction we take. In what 
follows we let M = {1, . . . m}. 

Definition 2.1. Let H — Uiez& where X C M. We introduce the new map 
fu ■ I — > I defined by 



We call the set H a hole and the function fu '■ I — > I the open dynamical system 
generated by the (closed) dynamical system / : I — >• I over H . 

The partition £ will remain a Markov partition of the open dynamical system 
fjj ■ I — > I but the dynamics of the original system / : I I will have been 
modified such that each point in H is now a fixed point. 

For n > let 



The set X n (fH) consists of those points that escape through the hole H at time n 
while Y n (fii) are those points that escape through H before time n + 1. 

For the moment suppose fi is a probability measure on J, i.e. //(/) = 1. Then 
/i(X™(/#)) can be treated as the probability that an orbit of / enters H for the 
first time at time n and n(Y n (fH)) the probability that an orbit of / enters H 



2. Open and Closed Dynamical Systems 




X n {f H ) 



{x G I : f\x) G H, f k (x) <£H, < k < n} 
{xel: f%{x) G H, f%{x) £H, < k < n}; and 



Y n {f H ) 



{xel: f k {x) G H, for some k, < k < n} 
{xel: fff(x) e H, for some k, < k < n}. 
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before time n + 1. In this regard, 

p n (f H ) = i - KY n (.M) 

represents the probability that a typical point of / does not fall into the hole H by 
time n. For this reason, the quantity P n {fn) is called the survival probability at 
time n of the dynamical system Jh for the measure /z. 

One of fundamental problems in the theory of open systems is determining or 
finding ways of approximating and [i(Y n (f H )) for finite n > 0. In 

the following section we give exact formulae for these quantities in the case where 
fu is a piecewise linear function with nonzero slope and fi is Lebesgue measure. 
In section [3] we remove the assumption that fjj is piecewise linear and present a 
method for estimating /i(X n (/jy)) and n(Y n (j h)) for functions that are piecewise 
nonlinear. 

3. Piecewise Linear Functions 

As a first step, we consider those open systems fjj : I — > I that are linear when 
restricted to the elements of the partition £. More formally, suppose H — \J ieX & 
for some IcM. Let C be the set of all open systems fjj : I — >• I such that 

\f' H (x)\ = Ci > for x € & and i ^ X 

where each Cj € K. The set C consists of all open systems that have a nonzero 
constant slope when restricted to any £i <^ H. 

To each fjj G C there is an associated matrix that can be used to compute 
fi(X n (f H )) and fi(Y n (f H )). To define this matrix let 

(1) 6i = ein/- 1 (^)forl<*,i<m. 

Definition 3.1. Let fn G C where H = [J ieX (,i for some X C M. The matrix 
Ah G R mxm given by 

, fl/'^r 1 for xee«^0, iiT, , 
(A H ) ij = < . l<i,J<m 

I otherwise 

is called the weighted transition matrix of fn- 

Associated with the open system fn '■ I — > I there is also a directed graph 
F = (V,Eh,uj) with vertices V and edges Eh- For V — {vi, . . . ,v m } we let 
denote the edge from vertex Vi to Vj. 

Definition 3.2. Let fn-I—^I where H = Uiexf* f° r some X C M. We define 
Tff = (V, Eh) to be the graph with 

(a) vertices V = {v%, . . . , v m } and; 

(b) edges E H = {ey : c/fo) C d(/(&)), i £ X}. 
The graph r# is called the transition graph of fjj. 

The vertex set V of Th represents the elements of the Markov partition £ = 
{£i}"=i an( i the edge set E H , the possible transitions between the elements of £. 
Hence, e,j G only if there is an x G £i ^ # such that fn(x) G i.e. it is 
possible to transition from & to £j. We note that as H = is a possible hole, the 
original (closed) system /:/—►/ has a well defined transition graph which we 
denote by T. 

Note that the matrix Ah and the graph T# do not carry the same information 
as Th only designates how orbits can transition between elements of £ whereas Ah 
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Zi=H 6 & U 

Figure 1. The transition graph Tg (right) of the open system 
fH'-I—>I (left) in example [I] 

additionally gives each of these transitions a weight. However, the graph Th gives 
us a way of visualizing how orbits escape from the system, which will be useful in 
the following sections. An open system and its transition graph are demonstrated 
in the following example. 



Example 1. Let the function f : I I be the tent map 

' 2x 0<x<l/2, 
2-2x \j2<x<\ 



with Markov partition £ = {(0, 1/4], (1/4, 1/2], (1/2, 3/4], (3/4, 1]}. Here we con- 
sider the hole H = (0, 1/4]. The open system fu ■ I — > I is shown in figure^ (left) 
with the graph Th (right). 

As H = £i we emphasis this in Th by drawing the vertex v% as an open circle, 
i.e. as a hole. We note that the only difference between the transition graph T of 
f : I —¥ I and Th is that there are no edges origination from V\ in Y h . In this 
sense a hole H is an absorbing state since nothing leaves H once it enters. 



Let 1 = [1, . 



. , 1] be the 1 x m vector of ones and e# the m x 1 vector given by 
I if i G X 



(ejy)< 







otherwise 



Theorem 3.3. If fu G C and n > then 

(2) M* b (/h)) = l^W; and 



(3) 



»(Y n (f H )) = l(J2A i H )e H . 



i=0 



For the matrix B G R mxm let a(B) and p(B) denote the set of all eigenvalues 
and spectral radius of B respectively. We then have the following corollaries to 
theorem 13.31 



Corollary 1. // fn G C and p(An) < 1 then 



KY n (.M) = HI - A H )-\I - A n H +1 )e H - and 



lim f i(Y n (f H )) = l(I-A H )- 1 e H . 
where I is the identity matrix. 
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Corollary 2. Suppose f H £ £. // < p(A H ) < 1 then lim P n (f H )=0. If 

n—too 

p(Ah) = i/ien p(Y n (fH)) = 1 /or some n < oo. 

A matrix B £ fl£ mxm j s called defective if it does not have an eigenbasis, i.e. if 
there are not enough linearly independent eigenvectors of B to form a basis of K. m . 
A matrix with an eigenbasis is called nondefective. 

Corollary 3. Let fn £ C and suppose the matrix Ajj is nondefective with eigen- 
pairs {(Ai, vi), . . . , (A^, v,t)} with no eigenvalue equal to 1. Then e# = J2i=i 
for some C\, . . . , Cfc £ C and 



(4) 
(5) 

where s.; = lv; 



m*"(/h))=X> s * a " 



i=l 



p(Y n {f H )) = Y J C i s i 



1 - A 



n+l 



1-A,: 



Example 2. Let the function f : I —> I be the tent map considered in example^ 
and let H=(0,l/4]- As fu S £■ one can calculate that fjj has the weighted transition 
matrix 


1/2 1/2 
1/2 1/2 
1/2 1/2 



.4 



H 



The matrix Ah is nondefective as its eigenvalues a (Ah) = { 
corresponding respectively to the linearly independent eigenvectors 



_ r i+75 l-VE 



0,0} 





■ o ■ 




" " 









" -1 " 


Vl = 


1+V5 
4 

1 + V5 


,v 2 = 


i-Vs 

4 

1-V5 


,v 3 = 




-1 


,v 4 = 


1 





4 




4 




1 









1 




1 







Since the vector en = [1/4,0,0,0] T can be written as 



V5 



= Vl 



20 + 20y/5 8y/5 
equations and ^ in corollary^ imply 

p(X n (f H )) = -1(5 + v / 5)A™ + -1(5 



-v 2 



rv 3 



r v 4, 



40 



^"(/h)) = 1:(5 + v / 5)- 



40 

v"+l 



40 

= 1 - 



1-Ax 



- V5)X%; and 
40 (5 ^ 5) T 



A^ +1 



1 

7! 



A" + - 



1 



\n+l 
A 2 ■ 



7Vo£e £/ia£ as p(.A#) < 1 then Urtin-^oo P n (fn) = 0. Hence, the probability of 
surviving indefinitely in this system for a typical x £ I is in fact zero. This can be 
seen in figure^ where both p(X n (fu)) and p(Y n (f}j)) are plotted. 
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4. Nonlinear Estimates 

We now consider the open systems fjj ; I — > I where / is allowed to be a 
nonlinear but differentiable function when restricted to the elements of £. The 
formulae we derive in this section allow us to give upper and lower bounds on 
ti{X n {f H )) and n{Y n {f H )) for any finite time n > 0. 

Suppose H = Uigi £» f° r some X C M. Let Af be the open systems /#:/—>•/ 
such that 

inf \f' H {x)\ > for and 
sup < oo for 

To each fjj € Af there are two associated matrices similar to the weighted transition 
matrix Ah defined for each open system in C 

Definition 4.1. Let fn G Af where H = Uiei^* ^ or some ^ C M. The matrix 
A H el"™ is defined by 

inf l/'Oz)!" 1 for te^m, iil, 
( ••'!;/)-, <( ' ^ l<*,J<m. 
otherwise 

Similarly, define the matrix A H € R mxm by 

f sup l/'^r 1 for fo^fl, i?Z 
(A H )ij = < x ^i l<i,j<m. 
I otherwise 

For / H ejV and n > let 

X n (f H ) = lA n H e H and X n (f H ) = lA H e H ; 

m m 

Y n (f H ) = 1(J2£h)ch and F"(.fe) = l(^^ ff )e H . 

Theorem 4.2. If f H e Af and n > then X n (f H ) < fi(X n (f H )) < and 
Z"(/if)<M(y n (//f))<F"(/if)- 
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Theorem 4.2 allows us to bound the amount of phase space that escapes through 
H at time n and before time n + 1. If the matrices A H and Ah are nondefective 
then we have the following result similar to corollary [3] 

Corollary 4. Let fn G J\f and suppose both A H and Ah are nondefective with 
eigenpairs {(Aj, VjJ, . . . , (A fc , v fe )} and {(Ai, Vi), . . . , (Afe, v*;)} luiift no eigenvalue 
equal to 1 . Then for each n > 



(0) 
(7) 



1 - A 



n+l 



1 - A, 



i=i 



< 



A" +1N 



l-Ai 



where s { = ly^, Sj = lv,;. e# = c i v i and e# = CiVj 



i=i 



The upper and lower bounds given in ([6|) are and X (/#) respectively. 

The upper and lower bounds given in (M) are YJ l (fu) and Y (fu) respectively. 



< x< 1/2 



Example 3. Consider the function g : I —± I given by 
_ j ^x - 2lx 2 + 28x 3 

~\t0-- x )- 21 ( X - x ) 2 + 28 ( X - x ) 3 V2 < x < 1 

with Markov partition £ = {(0, 1/4], (1/4, 1/2], (1/2, 3/4], (3/4, 1]} and H = (0,1/4]. 
TTie function g : I —> I can be considered to be a nonlinear version of the tent map 
f : I I in example^ In fact both systems have the same transition graph (see 
figure^. 

For g H : I -> I one can compute £23 = (V 4 ;- 44 ]> £23 = (-44,1/2], £23 = 
(1/2, .55], 63 = (-55, .3/4], £23 = (3/4, .94], and £23 = (.94, 1]. From this we find 





0.29 2/11 

0.29 2/11 

2/11 0.29 



and Ah = 






0.29 





4 0.29 

4 0.29 

4 
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n— axis 



Figure 4. The upper bounds X {g H ) and lower bounds 

for /i(X n (gH)) are shown for the open system gn ■ I — > I in 

example [3] 



As A H and Ah are nondefective then using corollary^ one can compute that 

(8) 0.15A™ + 0.15A2 - 0.6 < X n (g H ) < -0.02A™ + 0.02A^ + .25 

where A x = 0.41, A 2 = —0.12, Ai = 4.27. and A2 = —0.27. Plotting the inequalities 
in yields the picture in figure fJJ Here the shaded area indicates the region in 
which (i(X n (gH)) must lie. 



5. Improved Escape Estimates 

In this section we define a delayed first return map of an open system fn E A/", 
which we will use to improve the escape estimates given in theorem |4.2| A key step 
in this procedure is to choose a particular vertex set of Th over which this map 
will be defined. This requires that we know the cycle structure of Th- 

A path V in the graph Th = {V,Eh) is an ordered sequence of distinct vertices 
V\, . . . , Vk £ V such that e^+i € E for 1 < i < k — 1. If the vertices v\ and Vk are 
the same then V is a cycle. If S C V where V is the vertex set of the graph we will 
write S = V\S. 

Definition 5.1. Let H = U ie x& for some ^ C M and let T H = (V,E H ). The set 
S C V = {vi, . . . , v m } is an open structural set of Th if V, G S for i£l and Th \s 
has no cycles. 

Structural sets were first defined in [B] where they were used to gain improved 
estimates of a dynamical network's stability. Later in [7] they were used to improve 
the eigenvalue estimates of Gershgorin et. al. Here, our goal is to extend their use 
to improve our estimates of [i(X n (f H )) and /j,(Y n (f H )). 

For the open system fn :/—>•/ we let st(Tn ) denote the set of all open structural 
sets of r#. If S G st(r^f) we let 1$ = {i E M : vi E S} be the index set of S and 
£s = U e x s &- 

Definition 5.2. Let S £ si(r^). For x E I wc let 7(2;) = i ?i ■ ■ - it where ij = k 
if fn(x) E £fc and t is the smallest k > such that f%{x) E £5. The set 

= {7 : 7 = 7(a0 for some x E I \ H} 

are the admissible sequences of with respect to S. 
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For x € I we say 7(3;) = ioi\ ■ ■ - it has length |7(x)| = t. The reason |7(s)| < 00 
is that the graph Tjj\g has no cycles. Hence, after a finite number of steps f H {x) 
must enter £5. 

Definition 5.3. Let S G st(Tff). For x E I and k > we inductively define 
x k+ i = Rf s {x k , ...,x ) where 

f l H {xk)l (xk) if Xk-i = x k for each < i < \j(x k )\ - 1, 
Xfc otherwise. 



Xk+l 



The function : I k+1 — > / is called the delayed first return map of fu with 
respect to 5. The sequence Xq,Xi,X2) ... is the orbit of Xo under -R/s. 

If T = max. x ei \ j(x)\ then strictly speaking x k +i = Rfs(x k , ■ ■ ■ ,x T ) for some 
t < T. The map Rfg acts almost like a first return map of fjj to the set £g. 
The difference is that a return to £s does not happen instantaneously (as it would 
happen in the case of a first return map) but is delayed so that the trajectory of a 
point under fu and Rfs coincide after a return to £g. 

For n > we let Rf§(xo) = x n and define 

X n (Rfs) = {xel: Rfij(x) G JT, £ JT, < k < n}; and 

Y n (Rfs) = {x G I : fl/JO) G ff, for some k, < fc < n}. 
Lemma 5.4. If S £ st^(Tf H ) and n > ^ erl = X n (Rf s ). 

Proof. For x £ I let 7(2:0) = «o*i ■■• where ij = k if fjf{xo) G Choosing 
5 G si(r#) let 7 S (a; ) = where = k if Rf s (x ) G £*• Let t > be 

the smallest number such that it G I5. Then 7(2:0) = *o*i ■ • ■ it & n d definition 5.3 
implies Rf s (xo) — f t H (xa). Therefore, i t = it where t G Is- 

Continuing in this manner it follows that ij = £j for each j elg. Since H C £g 
the point Xo, if it escapes, will escape for both fu and Rfs at exactly the same 
time. This implies the result. □ 

The major idea in this section is that one can use Rfs to study the escape of 
fn through H. However, the weighted transition matrix of Rfs ■ I k — > I cannot 
be defined in the same way that we have defined either A H or Ah- To define a 
transition matrix for Rfs we require the following. 

For S G st(T H ) let 

M s = M U {7; % : 7 G Cl s , < i < |7|} 

If 7 = io . . . , i t we identify the index 7; with i and the index 7; t with i t . We 
also let £ 7 = {x G I : j(x) = 7} for each admissible sequence 7 G f2s, which simply 
extends our notation given by ([I]) in section [3j 

Definition 5.5. For S G s£(r#) let A s be the matrix with rows and columns 

indexed by elements of Ms where 

(9) 

{inf |(/ l7l (a;))'| _1 iN = 75 111 - h 3 = 75 \~f\, for some 7 G fls 
1 if i = 7; fe — 1, j = 7; fc, A; 7^ |7|, for some 7 G Sis 

otherwise. 

We call A s the lower transition matrix of Rfs- The matrix As defined by replacing 
the infimum in ([£]) by a supremum is the upper transition matrix of Rfs ■ 
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Let Is be the 1 x \Mg\ vector given by 

' 1 if i e M, 



(Is), = 



otherwise. 



Let eg be the \Ms\ x 1 vector given by 

_ J/xfe) if iGl, 



( e s); = 



|0 otherwise. 



Lastly, for n > let 

X"(i?/ S ) = l s A^es and X"(i?/ 5 ) = l S ^e s ; 

n n 

T"(i?/s) = 1 S (]T A l s )e s and = ls(£4)es- 

i=0 i=0 

Using these quantities we give the following improved escape estimates. 

Theorem 5.6. Let Jh E N and suppose S E st(Tji)- lfn>0 then 

X n (f H ) < X n (Rf s ) < KX n {f H )) < X n (Rf s ) < X n (f H ); and 
Y n {f H ) < Y n (Rfs) < KY n (f H )) < Y n (Rfs) < Y n (f H ). 



Theorem |5.6| together will lemma [574] imply that the escape of fjj through H is 
better approximated by considering any of its delayed first return maps Rfs than 
fjj itself. We now give a proof of theorem |5.6| 



Proof. For S G si(T#) suppose i e M\I and j e I. Then 

(A We 1 _ /inf^- l/'WI-Vfe) if &^0, < , . f 
I otherwise 

To show that a similar formula holds for larger powers of A s suppose k G Is- If 
ik, kj G ilg then 

(A s ) 4fc (A 5 ) fc3 (e s ), = inf If'ix)]- 1 inf |/'(z)rVfe) 

<^{x G 6 : /ir(aO G ZkJ 2 H ( x ) e £?}• 
If either ife ^ fi s or fcj ^ £7 S then (As)t*Gis)jy ( e s)j = °- 



Suppose fc G M \ I5. If ifcj G fig then ikj; 1 G M5 and 

=1- inf K/ 2 (as))'| — 1 / 

<M{z G 6 : /it (a) € fk, /bW € 



(4s)i,«JiiUs)ifcJ ; ij(es)i=l- mf |(/V))'rV&) 



If ifej ^ ils then ij'fc; 1 ^ Ms- Since 

G4|)ij(e S )j = J] (£s)ik(A s )kjKQ 

k£M s 

= {As)ik(As)kjv(€j)+ (As)i,ikr,i(As)iky,i,j[J>(€j) 

< Y V{x € & : /ff(a;) G a, /h 0*0 € 
feex s u(Af\i s ) 
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Continuing in this manner it follows that 
(10) (^)«(es)i 

for i € M \X, j G X, and n > 1. Since (es)j = if j ^ X then for n > 1 equation 
(10 1 implies 

lsA&s = E E (Mh(es) 3 < E ^ € 6 : G 
ieM jeM s ieM\x 

= VL{x€l\H:fi(x)€H}. 

As l s A^e s = fi(H) then X n (Rf s ) < fi(X n (f H )) < X n (Rf s ) for n > where the 
second inequality follows by using the same argument with the matrix As- 

To show that < X_ n (Rfs) we again suppose that i € M \ X and j e X. 

Then we have 

Cdff)tf(e*)j - | Q otherw . se - Ushteh 

For larger matrix powers we have 

E Uh)*^)^) = E in ; i/'^i -1 in / i/'wrVfe) 

— — ^Ck-i 

fcGM\I s feGAf\I s 

keM\x s 2:6 Ifc3 ifcjesis 
From this it follows that 

G^)y(e fl )j = E Uff)»-kUff)*i/*(6) + E (AhMAh)^^) 

kels k£M\X s 

- E (As)ik{As)kjfJ>(tj) + E Ui?)i.*fci;i(4ff)ifej;i,jM(6 - ) 
=(4s)«(es)r 

Again, continuing in this manner we have (A 7 H)iji e H)j < (^is)ij( e s)j f° r * G 
j e X, and n > 1. As 1# A^e# = = ls-A^eg then 

l H A n H e H = ^AhU^h), < E E (A n sU e s)i = isA^S 

ieM jel ieM jeM s 

for n > 0. Hence, < X n (Rf s ). 

By using the same argument with the matrix As we obtain the inequality 
X_ n (Rfs) < X_ n (fH)- The second set of inequalities in theorem 5.6 then follow, 
which completes the proof. □ 

Example 4. Consider the open system cjh : I I given in example [5| Observe 
that the vertex set S = {vi, 173,174} is an open structural set oJTh since v± is the 
vertex that corresponds to H and the graph Th\§ = { v 2} has no cycles (see figure 
i- 

The delayed first return map Rgs can be written as 

9H(Xk) ifxk £ $42, 

Rgs{xk) = { g 2 H {xk) if x k ,x k -i G $42, 
Xk otherwise. 
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£1 = H £ 2 £3 £42 £41 

Figure 5. The delayed first return map Rgs ■ I — > I in example 
[4] where Rgs is delayed on the set £42 shown in red. 



The map Rgs is shown in figure^ as a one- dimensional map but is colored red on 
£42 to indicate that in fact the system is delayed on this set. Specifically, gnix) is 
shown in blue and gjj (x) is shown in red where the trajectories of Rgjj stay in £42 
for two time-steps before leaving. 

To compute the upper and lower transition matrices of Rfa note that the sys- 
tem's admissible sequences are given by £ls = {23,24,33,34,424,423,41} implying 



(11) 



M s = {1,2, 3, 4, 424; 1,423;!}. 



FromQs we compute that £424 = (3/4, .85] and £423 = (-85, .94]. The other partition 
elements have been computed in example^ Using the order given in (11) we obtain 



A s = 



























.29 


2/11 














.29 


2/11 








2/11 











1 


1 











.26 











.26 















As = 



























4 


.29 














4 


.29 








29 











1 


1 











.73 











1.18 















The vectors Is and es are given by 

l s = [1, 1, 1, 1, 0, 0] and e s = [1/4, 0, 0, 0, 0, 0] T . 

Figure shows that X_ n (gH) < Xj L (Rgs) and X (Rgs) < X (gu) for a number of 
n— values and indicates the extent to which using the delayed first return map Rgn 
improves our estimates of pi(X n (g #)) . The shaded regions in these graphs represent 
the difference between these upper and lower estimates respectively. 



6. Conclusion 

Our results demonstrate that the theory of isospectral transformations can be 
effectively applied to open dynamical systems to obtain sharper estimates of a 
system's survival probabilities. Previous applications of this theory have allowed for 
improvements in the classical estimates of matrix spectra and in obtaining stronger 
sufficient conditions for the global stability of dynamical networks. We do not doubt 
that this theory can be applied to another problems as well. 

Concerning the results of the present paper it is easy to see that the improved 
estimates we obtained can also be found for a much broader class of open systems 



IMPROVED ESTIMATES OF SURVIVAL PROB. VIA ISOSPECTRAL TRANSFORMATIONS^ 




lower bounds upper bounds 

Figure 6. Comparison between X n (gH) and X" (Rqh) (left) and 
X (gn) and X (Rgn) (right) from examples pi and pi 



via the same technique. For instance, one could extend these techniques to higher 
dimensional systems. 
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